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ABSTRACT 


An  apparatus  for  the  direct  determination  of  the  second  inter¬ 
action  virial  coefficients  of  mixtures  has  been  constructedr  It  is  a 
modification  of  the  original  design  of  Knobler, 

Measurements  have  been  made  on  the  systems  nitrogen  -  carbon 
dioxide,  methane  -  carbon  dioxide,  and  methane  -  nitrogen*  The  first 
two  systems  were  studied  over  the  temperature  range  -60^0  to  120°C, 
while  the  latter  was  studied  over  the  temperature  range  -100°C  to  120°C, 
The  experimental  results  were  interpreted  to  yield  values  of 
This  requires  a  knowledge  of  the  virial  coefficients  of  the  pure 
materials*  For  methane  and  nitrogen  the  literature  values  of  the 
second  virial  coefficients  are  in  good  agreement*  For  carbon  dioxide, 
at  temperatures  below  0°C,  the  published  values  of  the  second  virial 
coefficient  are  in.  disagreement  by  as  much  as  15  percent.  Using  the 

present  apparatus,  values  of  the  second  virial  coefficients  of  pure 

o  o 

carbon  dioxide  were  determined  over  the  temperature  range  0  C  to  -60  C„ 
The  experimental  values  of  the  interaction  virial  coefficients 
B^2  ate  in  good  agreement  with  the  few  literature  data  available. 

The  results  were  compared  with  predictions  based  on  the  theory 
of  corresponding  states.  Both  the  correlation  of  Guggenheim  and 
McGlashan,  and  the  equation  of  Pitzer  and  Curl  were  used*  For  the 
methane  -  nitrogen  system,  the  Pitzer  and  Curl  correlation  gives  better 
agreement  with  this  work  than  the  correlation  of  Guggenheim  and  McGlashan, 
For  the  methane  -  carbon  dioxide  system,  the  Guggenheim  and  McGlashan 
correlation  gives  better  agreement  at  temperatures  below  300°K,  while 
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the  correlation  of  Pitzer  and  Curl  gives  better  agreement  at  tempera¬ 
tures  above  300°K.  For  the  nitrogen  -  carbon  dioxide  system,  the 
correlation  of  Guggenheim  and  McGlashan  is  in  better  agreement  with 
this  work. 

Two  common  molecular  potentials  were  considered,  the  Lennard- 
Jones  and  the  Kihara.  Published  values  of  the  parameters  for  the  pure 
components  were  combined  according  to  mixing  rules  and  the  interaction 
second  virial  coefficients  calculated  for  comparison  with  the  experi¬ 
mental  results.  For  the  methane  -  nitrogen  and  the  methane  -  carbon 
dioxide  systems,  the  Kihara  potential  gives  the  better  prediction. 

For  the  nitrogen  -  carbon  dioxide  system,  the  Lennard-Jones  potential 
is  in  good  agreement  with  this  work. 
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Chapter  I 


INTRODUCTION  AND  GENERAL  DISCUSSION 


Thermodynamic  data  are  directly  useful  in  process  design  and 

calculations.  They  are  also  valuable  for  testing  correlations  and 

methods  of  prediction.  Although  a  large  body  of  such  data  have  become 

available,  there  is  always  a  demand  for  reliable  data  on  new  systems. 

Such  data  may  also  be  used  to  calculate  other  thermodynamic  properties. 

Gas  phase  non-ideality  has  been  considered  from  a  number  of 

42 

points  of  view.  The  virial  equation  of  state  of  Onnes  ,  because  of 
its  firm  theoretical  basis,  has  often  been  used  for  calculations 
involving  non-ideal  gases.  It  is  of  the  form: 


RT[1  +  B(T)/V  +  C (T )  ,/V 2  + 


] 


(1) 


PV 


where  P  is  the  pressure,  V  is  the  molar  volume,  T  is  the  absolute  temp¬ 
erature,  R  is  the  gas  constant  and  B  and  C  are  virial  coefficients. 

The  importance  of  this  equation  is  due  to  the  fact  that  it  converges 
rapidly  up  to  moderately  high  densities.  It  is  generally  possible  to 
truncate  this  series  at  the  linear  term  at  low  densities.  For  mixtures 
of  gases  the  second  virial  coefficient  is  given  by: 


(2) 


similarly,  the  third  virial  coefficient  is  written  as: 


i  j  k  ijk 


(3) 
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A  large  number  of  B  values  have  been  reported  in  the  litera- 

,  14,  32  J  r  u  . 

ture  tor  pure  gases  and  measurements  of  the  properties  of  a 

number  of  large  molecules  have  been  performed  but  only  a  limited  number 
of  interaction  virial  coefficients  have  been  determined.  There  is  a 
clear  need  therefore  to  measure  the  interaction  second  virial  coef¬ 
ficients  over  a  wide  range  of  temperature  to  provide  the  data  for 
calculation  of  thermodynamic  properties  of  mixtures.  The  relation  of 
second  virial  coefficients  to  the  intermolecular  potential  4>  (r) 
function  is  as  follows  for  a  spherically  symmetrical  molecule 


B  = 


*00 

[1  -  exp  (  -  ^)]db  <r) 


(4) 


where  r  is  the  distance  between  molecular  centers  and  b  =  2iTNr  /3. 

o 

It  can  be  clearly  understood  that  the  second  virial  coefficient, 
B,  is  dependent  on  temperature  on  the  basis  of  a  simple  physical 
picture.  At  low  temperatures  collisions  between  pairs  of  molecules 
are  strongly  influenced  by  the  long-range  attractive  intermolecular 
forces,  and  such  pairs  may  spend  considerable  time  in  the  vicinity  of 
one  another.  In  this  situation,  a  form  of  molecular  association 
exists  and  the  existence  of  these  transient  dimers  causes  a  reduction 
of  the  pressure  as  predicted  from  the  perfect  gas  law  so  B  is  negative. 
At  high  temperatures,  the  collisions  of  molecules  are  much  more  fre¬ 
quent  and  with  higher  energy;  also,  the  weak  attractive  forces  can  only 
slightly  influence  these  collisions.  Instead,  the  short-range  repul¬ 
sive  forces  of  the  molecules  dominate  in  energetic  collisions.  These 
effects  give  an  excluded  volume  and  produce  a  pressure  above  the 
perfect  gas  law  which  give  the  positive  value  of  B.  At  still  higher 
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temperature,  the  collisions  become  even  more  energetic  so  that  molecules 
feel  the  "softness"  of  each  other’s  cores,  and  the  excluded  volume 
effect  decreases.  This  is  the  reason  why  second  virial  coefficients 
decrease  again. 

It  is  well-known  that  in  the  simple  case,  the  naive  assumption, 
(B^  +  B,„)  means  that  there  is  no  pressure  change  on  mixing. 
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But  it  has  been  pointed  out  by  Guggenheim  that  the  naive  assumption 
is  in  contrast  with  experiment. 

Conventional  methods  for  the  study  of  the  equation  of  state  are 

difficult  and  required  great  care.  For  this  reason,  the  differential 

techniques  for  determining  virial  coefficients  have  been  developed. 

A  special  design  for  the  determination  of  interaction  virial  coeffi- 

28 

cients  has  been  developed  by  Knobler  .  In  this  method,  at  constant 
temperature,  two  identical  bulbs  are  filled  to  the  same  pressure  with 
different  gases.  A  Toepler  Pump  is  used  to  perform  the  mixing  and 
after  mixing,  the  pressure  change  is  read  on  an  oil  manometer  and  may 
be  directly  related  to  the  excess  second  virial  coefficient: 


E  =  B 


12 


2  (B11  +  B22) 


(5) 


This  method  has  been  applied  with  success  to  mixtures  of  the 

o  3  27  33 

permanent  gases  from  77  to  320  K  *  *  ,  The  restriction  of  this 

method  is  that  it  cannot  be  applied  to  mixtures  containing  condensable 
vapors.  Gases  cannot  be  studied  if  the  gases  are  appreciably  soluble 
in  oil  and  mercury.  Gases  with  the  pressure  above  the  saturation  at 
room  temperature  cannot  be  studied  because  the  mixing  of  gases  by 
Toepler  pump  has  to  be  at  room  temperature  to  reduce  the  vapor  pressure 


' 


9M 


4 


of  mercury.  The  method  used  in  this  work  is  a  modification  of 
28 

Knobler's-  design.  The  detailed  description  will  be  presented  in 
Chapter  II. 

Table  I  shows  a  review  of  all  the  directly  measured  experimental 
interaction  second  virial  coefficients  up  to  this  date. 

The  importance  of  the  virial  equation  of  state  lies  in  its 
theoretical  connection  with  the  forces  between  molecules.  These 
forces  can  be  considered  the  basic  unifying  element-  connecting  a  host 
of  widely  different  properties;  hence,  it  is  unnecessary  to  measure  all 
of  the  properties.  In  this  work  values  of  the  interaction  second 
virial  coefficients,  for  three  binary  pairs  have  been  determined 

over  a  wide  temperature  range.  This  makes  the  values  useful  for  cal¬ 
culations  involving  the  virial  equation  of  state  and  also  allows  the 
calculation  of  other  thermodynamic  properties  over  a  wide  range  of 


conditions . 
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TABLE  I 

Experimental  Interaction  Second  Virial  Coefficient  Data 


Year 

System 

Authors 

References 

1942 

N2  "  C02 

Edwards  and  Roseveare 

15 

1951 

N2  '  C02 

Lunbeck  and  Boerboom 

30 

1953 

13  binary  mixtures 

Michels  and  Boerboom 

40 

1953 

n2  -  co2 

Gorski  and  Miller 

17 

1955 

N2  -  C02 

Pfefferle  et  al. 

44 

1956 

c°2  -  He 

Cottrell  and  Hamilton 

9 

1956 

5  binary  mixtures 

Cottrell,  et  al. 

8 

1959 

1961 

13  binary  n 

N2  -  C02 

aixtures 

Knobler,  et  al. 

27 

1963 

n2  -  ch4 

CO.  -  CH,  , 
2  4  J 

He  -  N2 

Mason  and  Eakin 

31 

He  -  Ar 

Ar  -  N2 

» 

Magasanik 

33 

196  7 

55  binary  mixtures 

Brewer 

2 

1968 

15  Normal 
alkane  (C^ 
mixtures 

'  V 

Dantzler  et  aX„ 

10 

1969 

11  binary  mixtures 

Dantzler  and  Knobler 

11 

1969 

9  binary  mixtures 

Brewer  and  Vaughn 

4 

• 
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Chapter  II 

EXPERIMENTAL  METHODS 


In  this  chapter,  the  construction  of  the  equipment,  the  control¬ 
ling  and  measuring  instruments  used,  and  the  procedure  for  the  deter¬ 
mination  of  excess  second  virial  coefficients  are  described.  Deter¬ 
mination  of  second  coefficients  for  pure  materials  is  also  described. 

A.  General  Survey 

The  principle  upon  which  the  differential  method  of  the  excess 
second  virial  coefficient  is  based  is  the  fact  that  when  two  different 
gases  which  are  initially  at  the  same  temperature  and  pressure  are 
mixed,  a  pressure  change  will  occur.  Since  most  gaseous  mixtures  have 
negative  excess  second  virial  coefficients  in  this  experimental  tempera¬ 
ture  range,  the  pressure  always  increases  upon  mixing. 

It  has  been  pointed  out  in  the  previous  chapter  that  excess 
second  virial  coefficient  is  E  =  “  y(B.^  +  ®22^ '  ^he  Huantities 

determined  experimentally  are  the  filling  pressure  P,  and  the  pressure 
difference  after  mixing,  AP »  The  relation  of  AP  with  excess  second 
virial  coefficient  E  is  as  follows: 


E  =  2 


AP 

P 


(dl  ~  V 

4dld2 


RTdld2 


(dl  -  d2  )(Cmd1  c222d2)  3 

- 4dX - 4(dlGl  +  d2G2 


(6) 
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and  G2  are  excess  third  virial  coefficients  defined  by 
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G,  =  C.  .  .  -  t(2C. . .  +  C. . .) 
i  113  3  111  jjj' 


(7) 


Equation  (6)  is  reduced  upon  substitution  of: 


(  > 

p 

-  B 

(  \ 

P 

RT 

\  > 

ii 

RT 

+  (2B. .  -  C. .  .) 

11  111 


'P 13 


,RT, 


+  ... 


(8) 


for  the  density  of  the  two  pure  components.  This  then  gives, 


E  =  2RTAP  _  2 


(B11  "  B22)  /4 


E(B11  +  B22)/2 


+  3(G1  +  G2)/8 


,RT, 

V.  ✓ 


'  2  h(Bll  -  B22> 


^Clll  B11  ^  ”  (C222  "  B22  ) 


+ 

[e] 

*1 

(2BH  +  BnB22 


+  2B22  Clll  C222)  8 


G1(‘2B11  +  B22^ 


+  G2(B11  +  2B22} 


; 


(9) 


from  which  E  can  be  calculated  by  successive  computer  iterations. 
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Equation  (6)  was  originally  derived  by  Knobler  .  After  making 
substitutions,  Equation  (9)  has  a  different  form  from  Knobler's.  The 
detailed  derivation  is  shown  in  Appendix  D. 

B.  Experimental  Equipment  Construction 

The  apparatus  consists  of  six  main  parts,  namely, 

(1)  Three-compartment  stainless  steel  cylinders  with 


identical  volume. 
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(2)  A  temperature  bath  assembly  for  low  temperature  and 
high  temperature  operation. 

(3)  A  mercury  Toepler  pump. 

(4)  An  oil  and  two  mercury  manometers, 

(5)  A  precision  temperature  indicator. 

(6)  A  cathetometer . 

The  detailed  schematic  flow  diagram  of  the  apparatus  is  shown  in 
Figure  1.  Details  of  the  construction  are  as  follows. 

1.  Gas  containers 

The  three  stainless  steel  cylinders  were  made  of  three  2-1/2  in. 
inside  diameter  and  13  in.  long  pipes.  Each  of  the  compartments  has  a 
volume  of  about  1000  cc  and  the  difference  between  the  three  volumes  is 
within  0,1  percent.  These  three  cylinders  were  oriented  120  degrees  to 
a  center  small  tube  in  which  the  platinum  resistance  thermometer  is 
placed.  The  three  cylinders  were  connected  to  2  mm.  glass  capillary 
tubing  by  means  of  1/8  in.  O.D.  stainless  steel  tubing  outside  of  the 
bath  and  1/4  in.  O.D,  stainless  steel  tubing  inside  the  bath.  The 
stainless  steel  tubing  extended  nearly  to  the  bottom  of  each  cylinder. 
The  sections  of  the  tubing  within  the  cylinders  had  holes  drilled  in 
their  walls  to  make  the  flushing  and  charging  and  mixing  of  the  gases 
more  efficient. 

Two  of  the  valves  which  are  immersed  in  the  fluid  are  stainless 
steel  bellows  valves  (Hoke  Inc.,  type  4212  F4Y) .  The  other  valves  used 
are  screwed  bonnet  valves  (Whitey  Research  Tool  Co.). 

2.  Temperature  bath 

It  is  necessary  to  maintain  a  constant  temperature  within  0.01°C 
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1. 

Levelling  bulbs 

16. 

Sensor  for  the  refrigerator 

2. 

Mercury  manometer 

17. 

Refrigerator 

3. 

Stopcock 

18. 

Liquid  nitrogen  cooler 

4. 

Glass  tubing 

19. 

Nitrogen  vapor  exhaust 

5. 

Vacuum  pump 

control  valve 

6. 

Vacuum  release  valve 

20. 

Liquid  nitrogen  Dewar  (50  1 

7. 

Vacuum  valve 

21  . 

Pressure  reducer 

8. 

Constant  temperature  bath 

22. 

Air  pressure  line 

9  . 

Heater 

23. 

Air  release  line 

10. 

St  1 rrer 

24. 

Toepier  pump 

11. 

Stainless (steel  tubing  (1/8  in.  O.D.) 

25. 

Toepler  pump  stopcock 

12. 

Bellows  valves 

26. 

Oil  manometer 

13. 

Platinum  resistance  thermometer 

27,28. 

Oil  manometer  stopcocks 

14. 

Stainless  steel  cylinders  (gas 

29. 

Vent 

containers ) 

30. 

A  gas  container 

15. 

Sensor  for  the  heat  Input  controller 

31  . 

B  gas  coni  liner 

Figure  1 


Schematic  Diagram  of  the  Apparatus 
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during  the  experiment.  For  this  purpose  a  steel  bath  was  constructed 
and  4  in.  thickness  of  fibreglass  is  used  for  insulation.  The  net 
capacity  of  the  bath  is  about  170  liters.  The  platinum  resistance 
thermometer  is  placed  at  the  center  part  of  the  bath  to  measure  the 
absolute  temperature.  The  temperature  controller  used  a  thermistor 
probe  as  a  sensor  and  heater  of  750  watts.  For  the  high  temperature 
operation,  air  is  flowing  in  the  1/2  in,  O.D,  copper  coil  which  is 
installed  in  the  bath.  For  temperatures  below  room  temperature 
operation,  liquid  nitrogen  is  used.  The  bath  was  well-stirred  by  a 
single-stage  impeller.  The  fluids  used  are  Esso-10-W  oil  for  higher 
temperature  and  ethanol  for  the  temperatures  below  20°C. 

C .  Measuring  and  Controlling  Instruments 

The  instruments  which  are  used  to  measure  and  control  various 
valves  to  carry  out  this  experiment  are  described  briefly, 

1,  The  absolute  temperature  in  the  bath  was  determined  by  two 

platinum  resistance  thermometers  and  a  high  precision  inductive  ratio 

arm  bridge  (Automatic  Systems  Laboratories).  One  of  the  platinum 

resistance  thermometers  is  placed  in  the  bath  and  the  other  one  is 

placed  in  the  ice  dewar  as  a  reference.  The  indicated  ratio  from  the 

R1 

bridge  is  n  =  „ — — — —  ,  where  R  is  the  reference  platinum  resistance 

R  1  I  R  ry  -1- 

thermometer  and  R£  is  the  one  in  the  bath.  The  calibrations  of  the 
platinum  resistance  thermometers  are  presented  in  Appendix  A, 

2.  The  filling  pressure  is  measured  by  means  of  a  mercury  manometer. 
The  pressure  difference  after  mixing  is  measured  with  an  oil  manometer 
by  means  of  a  precision  cathetometer  (W.G.  Pye  &  Co.  Ltd.).  The  deter¬ 
mination  of  the  oil  density  is  presented  in  Appendix  A, 


11 


o 

3.  The  bath  temperatures  were  controlled  within  ±0.01  C  by  a  temp¬ 
erature  controller  (Hallikainen  Instruments) . 

4,  The  refrigeration  unit  (Cincinnati  Sub-Zero  products)  was  used 
to  maintain  constant  temperatures  overnight. 

D .  Experimental  Procedure 

When  the  temperature  equilibrium  between  the  bath  and  the 
cylinders  was  being  attained,  several  successive  flushings  of  the 
Toepler  pump  with  a  50-50  mixture  of  the  gases  at  2  atm.  under  study 
were  made.  The  stopcock  to  the  Toepler  pump  was  closed.  After  four 
successive  evacuations  and  flushings  with  one  of  the  components ,  two 
of  the  cylinders  were  charged  to  about  2  atm.  Then  with  three  more 
evacuations  and  flushings  of  the  third  cylinder,  it  was  charged  with 
the  second  component.  The  pressure  of  the  gases  in  the  cylinders  can 
be  read  from  the  mercury  manometer  and  approximately  equalized  by 
using  the  mercury  levelling  bulbs  with  the  aid  of  the  oil  manometers. 

After  one  hour,  the  gases  in  the  cylinders  attained  the 
equilibrium  temperature  with  the  bath.  The  stopcocks  of  the  oil  mano¬ 
meter  are  opened,  once  again  the  pressure  of  two  different  gases  can 
be  equalized  by  the  levelling  bulbs.  The  oil  manometer  is  isolated 
again,  then  the  stopcock  to  the  Toepler  pump  is  opened  and  the  mercury 
level  is  pressed  back  by  the  outside  air  line  to  the  fiducial  mark. 

The  small  differential  pressure  was  read  by  the  cathetometer  and  the 
reading  recorded  as  the  initial  differential  pressure. 

One  of  the  two  cylinders  which  had  the  same  gas  was  then  iso¬ 
lated  for  use  as  a  reference  pressure  and  the  oil  manometer  was 
isolated  again.  The  mixing  of  the  two  pure  gases  was  accomplished  by 
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operating  the  Toepler  pump  for  60  cycles.  The  operation  of  the  Toepler 
pump  was  performed  automatically  by  a  relay  unit  and  a  solenoid  valve 
on  the  air  line.  The  mercury  level  was  then  brought  back  to  the 
fiducial  mark  manually  and  the  Toepler  pump  stopcock  was  closed. 

After  opening  the  oil  manometer  stopcock,  the  differential  pressure 
after  mixing  was  read. 

To  make  sure  of  the  complete  mixing,  the  Toepler  pump  was  oper¬ 
ated  for  another  five  cycles  and  it  was  found  that  the  same  pressure 
difference  was  read  for  most  of  the  cases.  Then  the  reading  was 
accepted  as  the  final  differential  pressure  after  mixing.  The  actual 
differential  pressure  is  simply  the  difference  between  the  final  and 
initial  reading.  The  experimental  run  was  replicated  with  every 
temperature , 

A  small  part  of  the  gases  under  study  are  exposed  to  room  temp¬ 
erature.  Serious  error  will  occur  when  the  volume  of  gas  at  the  con¬ 
trolled  temperature  and  the  volume  of  gas  at  room  temperature  are  com¬ 
parable,  For  this  design,  only  0,47  percent  of  the  volume  of  the  gas 
is  exposed  to  room  temperature.  The  tubing  from  the  bath  to  the 
Toepler  pump  and  the  oil  manometer  was  made  as  short  as  possible. 
Stainless  steel  tubing  of  1/8  in,  0,D.  and  2  mm,  I.D.  capillary  glass 
tubing  were  used  for  this  purpose. 

This  experiment  cannot  be  carried  out  at  higher  pressures 
because  the  third  virial  coefficient  of  the  pure  component  will  affect 
the  differential  pressure  after  mixing  and  increase  the  experimental 


error . 
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E .  Experimental  Determination  of  the  Second  Virial  Coefficient  of 

Carbon  Dioxide 

The  second  virial  coefficient  of  carbon  dioxide  below  0°C  has 
been  determined  in  this  work  because  of  the  disagreement  of  the  data 
in  the  literature „ 

The  same  apparatus  which  has  been  used  for  measuring  interaction 
second  virial  coefficients  was  used,  It  was  accomplished  by  measuring 
the  change  in  pressure  that  occurs  between  helium  and  carbon  dioxide 
when  the  temperature  is  varied.  This  determination  permits  the  calcu¬ 
lation  of  the  change  in  the  second  virial  coefficients  with  temperature. 
One  value  of  the  second  virial  coefficient  at  the  reference  temperature 

must  be  known.  In  this  work,  B  =  -150  cc/gm-mole  at  0°C  was  used. 

lu  2 

The  virial  equation  of  state  is  used  for  the  calculation  and  the 
method  of  calculation  is  presented  in  Appendix  E. 

These  new  determinations  of  the  second  virial  coefficient  for 
pure  carbon  dioxide  are  given  in  Table  V. 

F .  Purity  of  Gases 

The  sources  and  purity  of  the  gases  used  are  given  in  Table  II. 
The  gases  were  used  without  further  purification. 
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TABLE  II 


Sources  and  Purity  of  the  Gases  Used 


1.  Nitrogen  (high  purity) 


Supplier:  Consumer's  Welding  Co. 

Typical  Analytical  Data: 


Nitrogen 

99.993  mole 

Oxygen 

20 

ppm 

Carbon  dioxide 

5 

ppm 

Hydrogen 

5 

ppm 

Moisture 

15 

ppm 

Argon 

100 

ppm 

Maximum  total  impurities 

40 

ppm 

(non-inert) 

2.  Methane  (ultra  pure) 

Supplier:  Matheson  Gas  Products 


Actual  Analytical  Data: 

Methane 

Carbon  dioxide 

Oxygen 

Nitrogen 

Ethane 

Propane 


99,97% 

<  10  ppm 

9  ppm 
24  ppm 
35  ppm 

<  5  ppm 


3.  Carbon  Dioxide 


Supplier:  Canadian  Liquid  Air  Co.  Ltd. 

Typical  Analytical  Data: 

Carbon  dioxide 
Sulfur  compounds 
Carbon  monoxide 


99.95%  by  volume 
0  ppm 
0  ppm 


15 


TABLE  II  (continued) 


Oil 

Balance 


4.  Helium 

Supplier:  Canadian  Liquid  Air  Co.  Ltd. 

Purity:  99,995% 


5  ppm 
inerts 
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Chapter  III 


EXPERIMENTAL  RESULTS  AND  THEIR  COMPARISON 


This  chapter  covers  the  experimental  results  and  the  comparison 

of  this  work  with  other  experimental  work  in  the  literature. 

The  second  interaction  virial  coefficients  are  calculated  from 

the  second  and  third  virial  coefficients  of  the  pure  gases.  In 

section  A,  the  virial  coefficients  of  the  pure  gases  are  listed. 

In  section  B,  experimental  results  of  the  excess  second  virial 

coefficient  and  interaction  second  virial  coefficients  for  three 

systems,  methane  -  nitrogen,  methane  -  carbon  dioxide,  and  nitrogen  - 

carbon  dioxide  are  reported.  The  comparison  with  other  experimental 

values  is  made  in  this  section. 

A.  Pure  Component  Virial  Coefficients 

1.  Second  and  third  virial  coefficients  of  methane 

Second  virial  coefficients  obtained  by  different  sets  of  workers 

are  all  in  good  agreement.  A  smooth  curve  through  all  the  data  was 

14 

prepared  by  Dymond  and  Smith  in  1969.  This  curve  agrees  well  with 

16 

that  reported  by  Goodwin  in  1970, 

The  reported  third  virial  coefficients  of  methane  are  not  in 
good  agreement  with  each  other.  For  temperatures  below  0°C,  very  few 

sets  of  the  data  can  be  compared.  The  values  used  in  this  work  are  the 

_  2,13,22,35,36,48 

smoothed  values  from  several  workers  in  the  literature 

The  listing  of  the  second  and  third  virial  coefficients  of  methane  used 
for  this  work  is  presented  in  Table  III. 


7.  < 
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Second 

TABLE  III 

and  Third  Virial  Coefficients 

of  Methane 

Temperature  (°K) 

B  (cc/gm-mole) 

2 

C  (cc  /gm-mole 

173.15 

-145.0 

4900 

193.15 

-116.0 

4300 

213.15 

-  94.5 

3860 

233.15 

-  78.0 

3450 

253.15 

-  65.0 

3050 

273.15 

-  54.5 

2700 

293.15 

-  45.0 

2600 

313.15 

-  37.5 

2390 

333.15 

-  31.0 

2250 

353.15 

-  26.0 

2140 

373.15 

-  21.0 

2120 

393.15 

-  16.5 

1880 
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2.  Second  and  third  virial  coefficients  of  nitrogen 

Values  of  the  second  virial  coefficient  of  nitrogen  which  have 

been  reported  generally  agree  to  within  the  estimated  experimental 

uncertainties.  A  smooth  curve  through  all  the  data  is  also  prepared 

14 

by  Dymond  and  Smith  .  The  third  virial  coefficients  obtained  by 
different  sets  of  workers  agree  well.  Smoothed  values  from  these 
workers^*  * ^  ^  are  used.  Table  IV  shows  the  second  and  third 

virial  coefficients  of  nitrogen  which  are  used  in  this  work. 

3.  Second  and  third  virial  coefficients  of  carbon  dioxide 

At  temperatures  above  200°K,  the  second  virial  coefficients 

5  12 

reported  by  Dadson  et  al.  *  are  significantly  less  negative  than 

7  34  39 

those  obtained  by  most  other  workers  *  *  .  It  is  not  possible  to 

say  which  set  of  data  is  the  more  accurate.  At  temperatures  below 

o  ••  49 

250  K  there  is  a  marked  discrepancy  among  the  results  of  Schafer  , 

7  2 

and  Cook  and  Brewer  .  The  values  of  the  interaction  second  virial 

coefficients  will  depend  on  the  pure  second  virial  coefficients  that 

have  been  chosen.  The  values  between  270  and  400°K  for  this  work  are 

39 

the  smooth  curve  based  on  Michels  and  Michels  ,  MacCormack  and 
34  2 

Schneider  and  Brewer  because  these  data  are  in  good  agreement.  The 
values  below  270°K  had  been  determined  by  this  work  and  are  in  good 
agreement  with  Brewer.  Figure  2  shows  the  comparison  and  the  smooth 
curve  used  in  this  work. 

The  reported  third  virial  coefficients  of  carbon  dioxide  is 

even  more  scarce  especially  for  temperatures  below  0°C.  A  smooth  curve 

2  39  5 

from  Brewer  ,  Michels  and  Michels  and  Butcher  and  Dadson  is  used  in 


this  work. 
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TABLE  IV 


Second  and 

Third  Virial  Coefficients 

of  Nitrogen 

Temperature 

(°K) 

B  (cc/gm-mole) 

C  (cc^/gm-mo 

173.15 

-51.0 

2250 

193.15 

-39.0 

2000 

213.15 

-29.0 

1900 

233,15 

-21.5 

1800 

253.15 

-15.5 

1700 

273.15 

-10.5 

1500 

293.15 

-  5.8 

1250 

313.15 

-  2.0 

1250 

333.15 

+  1.0 

1250 

353.15 

+  3.8 

1200 

373.15 

+  6.0 

1200 

393.15 

+  8.2 

1150 

( cc  /  g m  -  mole  ) 
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TEMPERATURE  (°K) 

Figure  2.  Second  Virial  Coefficients  of  CO2 
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The  values  of  second  and  third  virial  coefficients  of  carbon 
dioxide  used  in  this  work  are  listed  in  Table  V. 

B.  Interaction  Second  Virial  Coefficients 

The  experimental  results  of  interaction  second  virial  coeffi¬ 
cients  of  the  three  systems  are  presented  in  this  section.  The  raw 
data  which  show  the  filling  pressure,  the  differential  pressure  after 
mixing,  the  excess  quantity,  E,  the  second  and  third  virial  coeffi¬ 
cients  of  the  pure  components ,  and  finally  the  interaction  second 
virial  coefficients  are  presented  in  Appendix  B. 

1.  Methane  -  Nitrogen  System 

The  excess  second  virial  coefficient  and  interaction  virial 
coefficients  of  methane  -  nitrogen  mixtures  at  individual  tempera¬ 
tures  are  given  in  Table  VI. 

The  experimental  results  from  other  workers  and  calculated 

interaction  virial  coefficients  from  the  extrapolation  of  P-v-T  data 

of  other  workers  are  shown  in  Table  VII.  Figure  3  shows  that  Brewer's 

experimental  result  and  the  results  derived  from  P-v-T  data  of  Keyes 
24 

and  Burks  are  in  good  agreement  with  this  work. 

No  comparisons  can  be  made  for  temperatures  below  0°C  because 
no  data  are  available  in  literature. 

2.  Methane  -  Carbon  Dioxide  System 

Table  VIII  shows  the  excess  quantities  and  the  interaction 

second  virial  coefficients  of  methane  -  carbon  dioxide  mixtures  of  this 

47 

work.  The  P-v-T  data  of  Reamer,  et  al.  have  been  extrapolated  to 

give  B,,  then  interaction  second  virial  coefficients  of  methane  - 
M 

carbon  dioxide  were  calculated  and  are  presented  in  Table  IX.  The 
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TABLE  V 

Second  and  Third  Virial  Coefficients  of  Carbon  Dioxide 


Temperature 


(°K) 


B  (cc/gm-mole) 


C  (cc  /gm-mole  ) 


213.15 

-262.0 

1000 

233.15 

-214.0 

4700 

253.15 

-180.0 

5600 

273.15 

-150.0 

5750 

293.15 

-127.0 

4600 

313.15 

-110.0 

4600 

333.15 

-  95.0 

4400 

353.15 

-  83.5 

3900 

373.15 

-  73.0 

3400 

393.15 

-  62.0 

2900 

t 


■ 

23 


Excess 

and  Interaction 

TABLE  VI 

Second  Virial  Coefficients  of  the 

Methane 

-  Nitrogen  System 

Temperature 

(°K)  E 

(cc/gm-mole) 

(cc/gm-mole) 

173.15 

10.67 

-87.32 

193.15 

9.45 

-68.05 

213.15 

7.84 

-53.91 

233.15 

7.38 

-42.36 

253.15 

6.83 

-33.42 

273.15 

5.78 

-26.72 

293.15 

5.81 

-19.58 

313.15 

4.62 

-15.12 

333.15 

4.41 

-10.59 

353.15 

3.95 

-  7.14 

373.15 

2.31 

-  5.14 

393.15 

1.68 

-  2.47 
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TABLE  VII 

Interaction  Second  Vlrial  Coefficients  of  the  Methane  - 
Nitrogen  System  From  the  Literature 


Temperature  (°K) 

2  31  2  A 

Brewer  Mason  and  Eakin  Keyes  and  Burks 

273.15 

-25.72  -26.59 

288.50 

-16.8 

313.15 

-16.11 

353.15 

-  5.66 

25 


TEMPERATURE  (°K  ) 


Figure  3.  Comparison  of  Experimental  Results 
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Excess 

and  Interaction 

TABLE  VIII 

Second  Virial 

Coefficients 

of  the 

Methane  - 

Carbon  Dioxide 

System 

Temperature 

(°K) 

E  (cc/gm-mole) 

B12 

(cc/gm-mo 

213.15 

41.51 

-136.74 

233.15 

31.05 

-114.94 

253.15 

25.95 

-  97.55 

273.15 

21.00 

-  81.24 

293.15 

18.74 

-  67.25 

313.15 

16.42 

-  57.33 

333.15 

14.51 

-  48.49 

353.15 

11.62 

-  43.14 

373.15 

10,43 

-  36.56 

393.15 

5.98 

-  33.26 
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interaction  second  virial  coefficients  determined  by  other  workers  in 

the  literature  are  also  shown  in  Table  IX. 

Zaalishvili”*’*'  presents  values  of  interaction  second  virial 

coefficients  for  methane  -  carbon  dioxide.  From  the  P-v-T  data  of 
47  51 

Reamer,  et  al.  ,  Zaalishvili  calculated  values  for  the  second  virial 

coefficient  of  methane  -  carbon  mixture.  The  interaction  second  virial 

coefficient  values  are  calculated  by  using  second  virial  coefficients 

of  methane  from  his  previous  work‘d ,  and  carbon  dioxide  from 

34 

MacCormack  and  Schneider  .  These  values  are  also  given  in  Table  IX. 

Figure  4  shows  the  comparison  of  the  values  presented  in  Table 
IX  with  this  work.  Interaction  second  virial  coefficients  derived  from 
Reamer,  et  al  by  this  work  are  in  good  agreement  with  this  experimental 
work.  Due  to  the  different  values  of  pure  second  virial  coefficients 
of  methane  used  by  Zaalishvili,  the  interaction  second  virial  coef¬ 
ficients  of  Zaalishvili  are  slightly  different  from  those  derived  by 
this  work.  They  are  in  good  agreement  with  this  experimental  work 
because  the  second  virial  coefficient  of  methane  used  is  little 
different  from  that  used  in  this  work. 

At  temperatures  below  0°C,  there  are  no  reported  results  from 
the  literature,  so  no  comparison  can  be  made. 

3.  Nitrogen  -  Carbon  Dioxide  System 

The  experimental  results  of  this  work  are  listed  in  Table  X. 

Table  XI  presents  the  experimental  results  determined  by  other  workers, 

20 

and  the  extrapolated  results  from  the  P-v-T  data  of  Haney  and  Bliss 
Figure  5  shows  that  results  derived  from  P-v-T  data  of  Haney  and  Bliss 
are  in  very  good  agreement  with  this  work. 
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TABLE  IX 

Interaction  Second  Virial  Coefficient  of  the  Methane  - 
Carbon  Dioxide  System  From  the  Literature 


Temperature  Reamer  et  al 


47  V  1  •  t,  -51 

Zaalishvili 


Brewer  Mason  and  Eakin 


31 


273.15 

288.50 

311.0 

-57.42 

-63.2 

344.0 

-43.39 

-48.2 

377.0 

-32.36 

-36 . 6 

411.0 

-22.08 

-19.8 

-77.55 


-59.40 


(  cc  /  gm  -  mole  ) 
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TEMPERATURE  (°K) 

Figure  4.  Comparison  of  Experimental  Results 
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Excess 

TABLE  X 

and  Interaction  Second  Virial 

Coefficient 

of  the 

Nitrogen  -  Carbon  Dioxide 

System 

Temperature 

(°K)  E  (cc/gm-mole) 

B12 

(cc/gm-mole) 

213.15 

50.35 

-95.14 

233.15 

33.75 

-83.99 

253.15 

28.81 

-68.93 

273.15 

22.96 

-57.28 

293.15 

20.70 

-45.70 

303.15 

20.03 

-40.98 

313.15 

18.74 

-37.26 

323.15 

17.71 

-33.28 

333.15 

16.52 

-30.47 

348.15 

15.32 

-26.18 

363.15 

13.40 

-22.91 

373.15 

12.58 

-20.71 

393.15 

11.34 

-15.56 

1 

Interaction  Second  Virial  Coefficients  of  the  Nitrogen 
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Figure  5.  Comparison  of  Experimental  Results 
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For  temperatures  lower  than  0°C,  only  two  data  points  from 

2 

Brewer  are  available,  which  are  less  negative  than  this  work. 
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CHAPTER  IV 

THEORETICAL  PREDICTIONS  OF  INTERACTION  SECOND 

VIRIAL  COEFFICIENTS 


The  calculation  of  interaction  second  virial  coefficients  by 
means  of  the  principle  of  corresponding  states  and  molecular  poten¬ 
tials  are  performed  in  this  chapter.  For  corresponding  states, 

interaction  second  virial  coefficients  calculated  from  the  Guggenheim 
19  45 

and  McGlashan  and  Pitzer  and  Curl  correlations  are  compared  with 
this  work.  For  molecular  potentials,  the  comparison  between  this  work 
and  values  calculated  from  the  Lennard-Jones  and  Kihara  potentials 
is  also  made. 

A.  Corresponding  States 

1.  Guggenheim  and  McGlashan  Correlation 

In  order  that  a  substance  conforms  to  the  principle  of  corres¬ 
ponding  states,  the  assumption  that  the  interaction  energy,  U, 
between  two  molecules  with  distance  apart  of  r  can  be  completely  de¬ 
fined  by  the  simple  relation  as 


U 

e 


r  'i 

r 


(10) 


Where  e  is  an  energy  and  a  is  a  length  both  characteristic  of  the  mole 
cule,  and  y  is  the  universal  function  for  all  the  molecules  concerned. 
Upon  substitution,  Equation  (4)  can  then  be  written  in  dimensionless 
form 


B 

2nNa3 


CO 

ey (r/a)' 

1  -  exp 

kT 

^  o 

l  J 

(ID 
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If  c  is  set  proportional  to  the  critical  volume  and  e/k  propor¬ 
tional  to  the  critical  temperature  T  ,  then 


B 

V 


(12) 


where  F  is  a  universal  function  of  the  reduced  temperature.  The  curve 

19 

for  equation  (12)  had  been  prepared  by  Guggenheim  and  McGlashan  from 
data  on  pure  material.  Equation  (12)  can  be  applied  to  a  binary  mix¬ 
ture,  with  semi-empirical  mixing  rules  for  the  characteristic  para¬ 
meters  as 


V 

c.  . 
3-3 


1 

8 


+  V 

c . 
J 


(13) 


and 


T 

c.  . 
ij 


Tc  > 

j 


1/2 


(14) 


The  interaction  second  virial  coefficient  is  then  calculated  by  means 
of  equation  (2) . 

2.  The  Correlation  of  Pitzer  and  Curl 

Pitzer  proposed  a  third  parameter  for  the  corresponding  states. 
He  defined  an  acentric  factor  go,  which  is  a  macroscopic  measure  of  the 
extent  to  which  the  force  field  around  a  molecule  deviates  from 
spherical  symmetry.  The  operational  definition  of  co  is  given  as 


10 


0.7 


co 


-log 


-  1.0000 


(15) 


. 
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g 

where  P  is  the  saturation  vapor  pressure.  When  applied  to  the  second 
virial  coefficient,  the  extended  theory  of  corresponding  states  asserts 
that  for  all  fluids  in  the  same  class, 


B 


to 


kT 


e 


(16) 


where  a  is  a  characteristic  molecule  size,  and  e/k  is  a  characteristic 
energy  expressed  in  units  of  temperature.  The  function  F^  depends  on 
the  acentric  factor  to, 

Pitzer  and  Curl  have  rewritten  equation  (16)  by  replacing  a 
and  e/k  with  macroscopic  parameters  as 


BP 

c 

RT 


=  F 


to 


T 

l  CJ 


(17) 


and  investigated  the  dependence  of  F  on  acentric  factor  at  constant 


to 


reduced  temperature.  Within  the  limits  of  experimental  information 


then  available,  F  is  a  linear  function  of  to,  thus: 


to 


BP 

( 

RT 


=  F 


(0) 


'-j 


+  ooF 


(1) 


(18) 


where 


(0) 


T 

l  CJ 


=  0.1445  -  0.330 


^  y-1 
T 


-  0.1385 


T 


-  0.0121 


fT  )“3 


,  c. 
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(1) 


=  0.073  +  0.46 


rT  '-1 


-  0.0073 


v  cj 


T 

-8 


-0.50 


(  \  ““2 
'T 


-  0.097 


(r 


-3 


This  correlation  was  extended  to  mixtures  with  the  following  equations: 

(19) 


Z  R  T 
C12  C12 


'12 


V 


'12 


Z 

C12 

The  mixing  rules  as 

W12 

and  equation  (13)  and  (14) . 

3.  Comparison  of  Results 

The  calculated  results  of  interaction  second  virial  coefficients 
from  corresponding  states  of  Guggenheim  and  McGlashan  and  the  calcula¬ 
tion  of  the  Pitzer  and  Curl  correlation  are  shown  in  Tables  XII,  XIII, 
and  XIV  and  the  values  are  plotted  on  Figures  6,  7  and  8  for  the 
methane  -  nitrogen,  methane  -  carbon  dioxide  and  nitrogen  -  carbon 
dioxide  systems  respectively. 

For  methane  -  nitrogen  systems,  the  predicted  results  from 
Pitzer  and  Curl  correlation  are  slightly  more  negative  for  temperatures 
lower  than  350°K  but  less  negative  at  temperatures  higher  than  350°K 
in  comparison  with  this  work.  The  Guggenheim  and  McGlashan  prediction 
is  always  more  negative  than  this  work. 


0.271  -  0.08o) 


12 


(20) 


2(“l  +  “2) 


(21) 


« 
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TABLE  XII 


Interaction  Second 

Virial  Coefficients 

of  Methane  -  Nitrogen 

Mixtures  From  This 

Work,  Guggenheim  and 

McGlashan  Correlation 

and 

Pitzer  and  Curl  Correlation 

Temperature 

This  Work 

Guggenheim 

Pitzer 

173.15 

-87.32 

-89.8 

-88.28 

193.15 

-68.05 

-69.0 

-70.06 

213.15 

-53.91 

-56.0 

-56.01 

233.15 

-42.36 

-45.4 

-44.86 

253.15 

-33.42 

-36.9 

-35.82 

273.15 

-26.72 

-30.2 

-28.36 

293.15 

-19.58 

-23.6 

-22.09 

313.15 

-15.12 

-18.9 

-16.75 

333.15 

-10.59 

-14.2 

-12.16 

353.15 

-  7.14 

-  9.45 

-  8.17 

373.15 

-  5.14 

-  6.63 

-  4.67 

393.15 

-  2.47 

-  3.78 

-  1.58 
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TEMPERATURE  (°K) 

Figure  6.  Comparison  of  With  Values  Calculated  From 

Corresponding  States 
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Interaction  Second 

TABLE  XIII 

Virial  Coefficients 

of  Methane  -  Carbon 

Dioxide  Mixtures  From 

This  Work,  Guggenheim  and  MgGlashan 

Correlation 

and 

Pitzer  and  Curl 

Correlation 

Temperature  (°K) 

This  Work 

Guggenheim 

Pitzer 

213.15 

-136.74 

-158.0 

-164.41 

233.15 

-114.94 

-120.0 

-133.58 

253.15 

-  97.55 

-  99.0 

-109.71 

273.15 

-  81.24 

-  86.5 

-  90.73 

293.15 

-  67.25 

-  74.0 

-  75.32 

313.15 

-  57.33 

-  63.4 

-  62.59 

333.15 

-  48.49 

-  56.6 

-  51.92 

353.15 

-  43.14 

-  48.0 

-  42.87 

373.15 

-  36.56 

-  43.2 

-  35.11 

393.15 

-  33.26 

-  36.5 

-  28.39 

(  cc  /  gm  -  mole  ) 
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Figure  7.  Comparison  of  B12  With  Values  Caleulated  from 
Corresponding  States 
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TABLE  XIV 


Interaction  Second 

.  Virial  Coefficient  of  Nitrogen  -  Carbon 

Dioxide  Mixtures 

From  This  Work, 

Guggenheim  and  McGlashan 

Correlation  and 

Pitzer  and  Curl 

Correlation  (cc/gm-mole) 

Temperature  (°K) 

This  Work 

Guggenheim 

Pitzer 

213.15 

-95.14 

-92.0 

-96.39 

233.15 

-83.99 

-73.5 

-76.49 

253.15 

-68.93 

-62.5 

-60.88 

273.15 

-57.28 

-53.4 

-48.34 

293.15 

-45.70 

-45.0 

-38.09 

303.15 

-40.98 

313.15 

-37.26 

-37.8 

-29.56 

323.15 

-33.28 

333.15 

-30.47 

-32.1 

-22.38 

348.15 

-26.18 

353.15 

-28.4 

-16.26 

363.15 

-22.91 

373.15 

-20.71 

-23.0 

-10.99 

393.15 

-15.56 

-18.4 

-  6.41 
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Figure  8. 
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For  methane  -  carbon  dioxide  systems,  the  Guggenheim  and 
McGlashan  correlation  is  always  slightly  more  negative  than  this  work 
in  the  temperature  range  of  this  work.  The  Pitzer  and  Curl  correla¬ 
tion  is  much  more  negative  than  this  work  at  low  temperatures , 
becomes  closer  as  the  temperature  increases,  and  after  350°K,  it  is 
less  negative. 

For  the  nitrogen  -  carbon  dioxide  system,  the  Pitzer  and  Curl 
correlation  is  always  less  negative  than  this  work.  Guggenheim  and 
McGlashan  is  less  negative  at  temperatures  lower  than  310°K  and 
becomes  more  negative  than  this  work  at  higher  temperatures. 

B.  Molecular  Potential 

It  is  well  known  that  two  molecules  repel  each  other  when  they 
come  close  together  and  attract  each  other  when  they  are  far  apart. 
The  force  of  interaction  F  between  two  spherical  non-polar  molecules 
is  a  function  of  the  intermolecular  separation  r.  It  is  more  con¬ 
venient  to  use  the  potential  energy  of  interaction  (r)  rather  than 
the  force  of  interaction  F(r).  These  two  functions  are  simply 
related : 


F(r) 


d(j> 

dr 


(22) 


f  00 


<Kr) 


F(r)  dr 


(23) 


'  r 

Theoretical  considerations  suggest  that  the  repulsive  potential 
should  be  an  exponential  function  of  intermolecular  separation,  but  it 
is  more  convenient^  to  represent  the  repulsive  potential  by  an  inverse- 


power  Law  of  type, 


' 


45 


<Kr) 


a 

n 

r 


(24) 


where  a  is  a  positive  constant  and  n  is  a  number  usually  taken  to  be 
between  8  and  6. 

In  order  to  take  both  repulsive  and  attractive  forces  between 
nonpolar  molecules  into  account,  it  is  customary  to  assume  that  the 
total  potential  energy  is  the  sum  of  the  two  separate  potentials : 


4>  (r)  =  c)>  (r)  +  <f>  (r) 

repulsive  attractive 


a  8 

n  m 
r  r 


(25) 


where  a,  8,  n  and  m  are  positive  constants  and  where  n  >  m.  This 

41 

equation  was  first  proposed  by  Mie  and  was  extensively  investigated 
by  Lennard-Jones . 

Equation  (25)  gives  the  potential  energy  of  two  molecules  as  a 

function  of  their  separation  and  it  is  apparent  that  at  some  distance 
* 

r  ,  this  energy  is  a  minimum;  this  minimum  energy  is  designated  by 
& 

4>(r)  •  By  algebraic  rearrangement,  Mie’s  potential  can  be  written  as 


<t>(r) 


1 

,  n ,  mN n-m 
e (n  /m  ) 

n-m 


(26) 


where  e  =  -<Kr)  and  a  is  the  intermolecular  distance  when  4>  (r)  =  0. 
Figure  9  shows  the  shape  of  the  potential  function. 

The  relation  of  second  virial  coefficients  to  the  inter- 
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Figure  9.  Potential  Function  With  Adjustable 
Parameter 


47 


molecular  potential  function,  c|)(r)  which  is  shown  in  equation  (4)  is  as 
follows : 


B 


‘0° 

(1  -  exp[-<j)(r)/kT])db  (r) 
o 


1.  Lennard-Jones  Potential 

The  Lennard-Jones  potential  is  the  form  when  Mie's  potential  is 
rewritten  with  n  =  12  and  m  =  6,  then 


<P  ( r ) 


(27) 


This  equation  gives  what  is  probably  the  best-known  two-parameter 
potential  for  small  nonpolar  molecules.  In  Lennard-Jones'  formula, 
the  repulsive  wall  is  not  vertical  but  has  a  finite  slope;  this  implies 
that  if  two  molecules  have  very  high  kinetic  energy,  they  may  be  able 
to  interpenetrate  to  separations  smaller  than  the  collision  diameter  a. 
Potential  functions  with  this  property  are  sometimes  called  "soft- 
sphere"  potentials. 

21 

Force  constants,  e/k,  a  and  bQ  are  given  and  numerical  results 

based  on  Lennard-Jones  potential  are  available  for  second  virial  coef- 
21 

ficients  .  It  can  be  extended  to  mixtures  by  means  of  simple  empirical 
mixing  rules  as 

°12  =  7(01  +  °2)  (28) 


and 


(29) 
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2.  Kihara  Potential 

According  to  Lennard-Jones *  potential,  molecules  consist  of 

point  centers  surrounded  by  "soft"  (i.e.  penetrable)  electron  clouds. 

An  alternate  picture  of  molecules  which  was  proposed  by  Kihara  is  to 

think  of  them  as  possessing  impenetrable  (hard)  cores  surrounded  by 

penetrable  electron  clouds.  In  crude  mechanical  terms,  Kihara* s 

model  (for  spherically  symmetric  molecules)  considers  a  molecule  to 

35  26 

be  a  hard  billiard  ball  with  a  foam-rubber  coat.  Kihara  *  writes 
a  potential  function  identical  to  that  of  Lennard-Jones,  except  that 
the  intermolecular  distance  is  taken  not  as  that  between  molecular 
centers  but  rather  as  the  minimum  distance  between  the  surfaces  of 
the  molecular  cores.  The  Kihara  potential  with  three  parameters  is: 


(r)  = 


d>  (r)  =  4e 


a-2a 


r-2a 


12 


a-2a 


r-2a 


for  r  <  2a 


for  r  >  2a 


(30) 


where  2a  is  molecular  core  diameter. 
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Prausnitz  and  Myers  have  shown  that  the  calculation  of  second 
virial  coefficients  by  means  of  the  Kihara  potential  agrees  with  some 
experimental  data.  The  parameters  for  Kihara  potential  have  been 
determined  by  Prausnitz  and  Myers  from  precision  P-v-T  data.  When 
applied  to  mixtures,  the  interaction  second  virial  coefficient  was 
calculated  from  the  following  equation  according  to  the  Kihara 
potential , 


' 
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B  =  |nNa3[a*3  +  3(26)  a*2  F^T*)  +  3(23)  a*  F2(T*) 


+  (23)  F3(T*)]  [1  +  a*]~3 


(31) 


^  O 

in  which  T  =  kT/e  and  a*  =  a 


a-2a 


The  mixing  rules  are: 


12 


*  *  *  9 

=  (T  •  T  ) 

K  1  2  } 


(32) 


and 

a12*  =  2(al  +  a2  }  (33) 

Tables  of  F^(T  ),  F2(T  )  and  F^ (T  )  are  given  by  Hirschf elder , 

21 

Curtiss  and  Bird 
3.  Comparison  of  Results 

The  calculated  results  from  Lennard-Jones  potential  and  the 
computer  results  of  Kihara  potential  for  methane  -  nitrogen,  methane  - 
carbon  dioxide  and  nitrogen  -  carbon  dioxide  systems  are  given  in 
Tables  XV,  XVI,  and  XVII  respectively.  These  data  are  plotted  in 
Figures  10,  11  and  12  for  methane  -  nitrogen,  methane  -  carbon  dioxide, 
and  nitrogen  -  carbon  dioxide  systems  respectively  for  comparison  with 
the  results  of  this  work. 

For  the  methane  -  nitrogen  system,  the  Kihara  potential  has  a 
better  comparison  than  the  Lennard-Jones  potential  with  this  work. 

But  the  Pitzer  and  Curl  correlation  is  the  best  among  the  four 
theoretical  predictions  to  describe  this  system  with  the  temperature 
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TABLE  XV 


Interaction  Second 

Virial  Coefficients 

of  Methane 

-  Nitrogen 

Mixtures 

From  This 

Work,  Lennard-Jones 

and  Kihara 

Potentials 

(cc/gm-mole) 

Temperature 

(°K) 

This  Work 

Lennard-Jones  Kihara 

173.15 

-87.32 

-92.43 

-91.02 

193.15 

-68.05 

-73.79 

-71.46 

213.15 

-53.91 

-59.20 

-56.93 

233.15 

-42.36 

-47.56 

-45.90 

253.15 

-33.42 

-38.03 

-36.33 

273.15 

-26.72 

-30.08 

-28.95 

293.15 

-19.58 

-23.31 

-22.74 

313.15 

-15.12 

-17.63 

-17.11 

333.15 

-10.59 

-12.66 

-12.45 

353.15 

-  7.14 

-  8.34 

-  8.53 

373.15 

-  5.14 

-  4.55 

-  5.13 

393.15 

-  2.47 

-  1.19 

-  2.08 
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150  200  250  300  350  400 


TEMPERATURE  (°K) 

Figure  10.  Comparison  of  %  2  With  Values  Calculated  From 
Molecular  Potentials 
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TABLE  XVI 

Interaction  Second  Virial  Coefficients  of  Methane  -  Carbon 

Dioxide  Mixtures  From  This  Work,  Lennard-Jones  and  Kihara 

Potentials 

(cc/gm-mole) 


Temperature  (°K) 

This  Work 

Lennard-Jones 

Kihara 

213.15 

-136.74 

-158.06 

-153.74 

233.15 

-114,94 

-133.00 

-128.06 

253.15 

-  97.55 

-112.94 

-106.40 

273.15 

-  81.24 

-  96.41 

-  89.39 

293.15 

-  67.25 

-  82.49 

-  76.25 

313.15 

-  57,33 

-  70.60 

-  65.47 

333.15 

-  48.49 

-  60.49 

-  55.23 

353.15 

-  43.14 

-  51.71 

-  48.03 

373.15 

-  36.56 

-  43.85 

-  40.34 

393.15 

-  33.26 

-  37.01 

-  33.75 
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TEMPERATURE  (°K) 


Figure  11.  Comparison  of  With  Values  Calculated  from 
Molecular  Potentials 
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TABLE  XVII 

Interaction  Second  Virial  Coefficients  of  Nitrogen  -  Carbon 

Dioxide  Mixtures  From  This  Work,  Lennard-Jones  and  Kihara 

Potentials 

(cc/gm-mole) 


Temperature  (°K) 

This  Work 

Lennard-Jones 

Kihara 

213.15 

-95.14 

-95.63 

-92.08 

233.15 

-83.99 

-78.95 

-74.91 

253.15 

-68.93 

-65.35 

-61.07 

273.15 

-57.28 

-54.07 

-50.56 

293.15 

-45.70 

-44.61 

-40.66 

303.15 

-40.98 

313.15 

-37.26 

-36.42 

-33.07 

323.15 

-33.28 

333.15 

-30.47 

-29.46 

-25.89 

348.15 

-26.18 

353.15 

-23.32 

-19.98 

363.15 

-22.91 

373.15 

-20.71 

-17.99 

-16.23 

393.15 

-15.56 

-13.23 

-11.34 
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TEMPERATURE  (°K) 

Figure  12.  Comparison  of  B.^  With  Values  Calculated  from 
Molecular  Potentials 
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range  of  this  work. 

For  methane  -  carbon  dioxide  systems,  the  Kihara  potential  has 
a  better  comparison  with  this  work. 

For  nitrogen  -  carbon  dioxide  systems,  the  Lennard-Jones 
potential  has  a  better  comparison  with  this  work. 
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SUMMARY  AND  CONCLUSIONS 

1.  The  modification  of  the  original  Knobler’s  design  for  direct 
measurement  of  interaction  second  virial  coefficients  has  been  made. 
The  experiment  was  carried  out  at  2  atm.  pressure  in  the  temperature 
range  between  -100  C  and  120°C.  The  second  virial  coefficient  of 
carbon  dioxide  was  determined  between  0°C  and  -60°C  by  the  same 
equipment . 

2.  A  study  has  been  done  for  the  systems  of  methane  -  nitrogen, 
methane  -  carbon  dioxide  and  nitrogen  -  carbon  dioxide.  Comparison 
was  made  with  other  results  in  literature  by  direct  measurement  of 
interaction  second  virial  coefficient. 

3.  The  P-v-T  data  of  mixtures  of  methane  -  nitrogen,  methane  - 
carbon  dioxide,  and  nitrogen  -  carbon  dioxide  have  been  extrapolated 
to  obtain  B^,  and  were  then  calculated  by  means  of  equation  (2). 
Three  sets  of  results  were  in  very  good  agreement  with  this  work. 

4.  The  comparison  of  interaction  second  virial  coefficients  of 
methane  -  nitrogen,  methane  -  carbon  dioxide  and  nitrogen  -  carbon 
dioxide  between  this  work  and  correlations  based  on  the  principle 
of  corresponding  states,  those  of  Guggenheim  and  McGlashan,  and 
Pitzer  and  Curl  had  been  made.  For  the  methane  -  nitrogen  system, 
the  Pitzer  and  Curl  correlation  gave  better  agreement  with  this  work. 
For  the  methane  -  carbon  dioxide  system,  the  Guggenheim  and  McGlashan 
correlation  gave  better  agreement  at  temperatures  below  300°K,  while 
the  correlation  of  Pitzer  and  Curl  was  in  good  comparison  at  tempera¬ 
tures  above  300°K.  For  the  nitrogen  -  carbon  dioxide  system,  the 
correlation  of  Guggenheim  and  McGlashan  was  in  better  agreement  with 


[  ; 
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this  work. 

5.  Interaction  second  virial  coefficients  calculated  from  Lennard- 
Jones  potential  and  Kihara  potential  have  been  studied  and  the  compari¬ 
son  with  this  work  was  made.  For  the  methane  -  nitrogen  and  the 
methane  -  carbon  dioxide  systems,  the  Kihara  potential  gave  the  better 
prediction.  For  the  nitrogen  -  carbon  dioxide  system,  the  Lennard- 
Jones  potential  was  in  good  agreement  with  this  work. 

6.  Comparisons  were  difficult  to  make  because  the  reported  second 
virial  coefficients  were  not  consistent  in  the  literature,  especially 
second  virial  coefficients  of  pure  carbon  dioxide.  Results  which  used 
the  same  pure  second  virial  coefficients  as  this  work  derived  from 
P-v-T  data  of  other  workers  were  in  good  agreement  with  this  work. 

The  lact  of  agreement  between  the  experimental  results  and  those  cal¬ 
culated  from  the  principle  of  corresponding  states  and  from  molecular 
potentials  may  be  due  to  the  mixing  rules  used  to  calculate  the  para¬ 
meters.  The  validity  of  these  mixing  rules  should  be  tested  by  cal¬ 
culating  the  mixture  parameters  from  the  present  experimental  data 
and  determining  the  form  of  the  mixing  rule. 
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APPENDIX  A 

CALIBRATIONS 


A-l 


I.  Calibration  Data  for  8163  model  platinum  resistance  Thermometers 

(Serial  Nos.  1724262,  1752229) 

This  thermometer  was  calibrated  for  use  in  the  range 
90.188  K  <  T  <  773 . 15°K  (-182.962°C  <  t  <  500.60°C)  on  the  Inter¬ 
national  practical  temperature  scale  of  1968. 

For  13 . 81°K  <  T  <  273.15°K 


W  =  Wcct  +  dw 

where  W  =  R/R(0°C)  and  Wcct  is  a  reference  function,  independent  of  any 
particular  thermometer. 

For  90 . 188°K  <  T  <  273.16°K 

dw  =  A4  t  +  C4  (t  -  100)  t2 
For  0°C  <  t  <  630 . 74°C 
t  =  tf  +  dt 


where 


and 


t'  is  determined  by 
w=l+At’+B  (tT)2 


dt*  =  0 .45  (t  ’  /100)  (• 


t' 

100 


-  DC 


tV _ 

419.58 


-  1)0 


tj _ 

630,74 


1) 


The  coefficient  in  the  above  relations  were  found  by  fixed-point  cali¬ 
brations,  i.e.  triple  point  of  water,  the  steam  or  tin  point  and  the 
zinc  point  using  continuous  thermometer  currents  of  1  and  2  mA. 

For  Serial  No.  1724262 


0mA  2  mA 


A  = 

3.98556 

X 

OO 

o 

1 — 1 

3.98553 

X 

OO 

o 

I — t 

B  = 

-5.87665 

X 

io"7 

-5.87775 

X 

o 

1 

A4  = 

7.2864 

X 

io"7 

6.7345 

X 

o 

r— I 

C4  = 

1.1035 

X 

io-14 

1.1382 

X 

io~ 14 

For  Serial  No 


1752229 


R(0°C) 


A- 2 


OmA  2  mA 


A  = 

3.98468 

X 

10  3 

3.98465 

X 

Kf3 

B  = 

-5.87540 

X 

p'- 

i 

o 

rH 

-5.87660 

X 

10“7 

A4  = 

-1.0425 

X 

10' 7 

-1.5988 

X 

io~7 

C4  = 

8.1559 

X 

10-15 

8.5156 

X 

10'15 

=  25.5  ohms  (approximately.) 


II.  Calibration  of  oil  density  for  manometer 


A- 3 


Kind  of  oil:  cencoil-B 

Temperature:  20°C 


Run  1 

Weight  of  bottle  (50  c.c.) 
Weight  of  bottle  and  oil 

weight  of  oil 

Run  2 

Weight  of  bottle 
Weight  of  bottle  and  oil 


20.6541  gm 
69.6642  gm 

49.0101  gm 

20.6584  gm 
69.6547  gm 


weight  of  oil 


48.9983  gm 


Average  weight  of  oil  (50  c.c.  at  20°C)  =  49.0042  gm 
Density  of  oil  =  0.9801  gm/c.c. 


APPENDIX  B 


EXPERIMENTAL  DATA 


PAW  DATA  F03  METHANE-NITROGEN  SYSTEM 
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APPENDIX  C 

ERROR  ANALYSIS 


C-l 


1.  The  temperatures  were  controlled  within  ±  0.01°C  for  each  temp¬ 
erature  level.  This  can  introduce  an  error  in  excess  second  virial 
coefficient  of  ±  0.06  percent  at  low  temperature  regions  and  ±  0.03 
percent  at  high  temperature  regions  for  methane  -  nitrogen  systems. 

For  methane  -  carbon  dioxide  system,  ±  0.4  percent  at  low  temperature 
region  and  ±  0.2  percent  at  high  temperature  region.  For  nitrogen  - 
carbon  dioxide  system,  ±  0.8  percent  at  low  temperature  region  and 

±  0.06  percent  at  high  temperature  region. 

2.  In  this  work,  pressure  reading  is  within  0.05  mm  oil  and  using 
first  term  of  equation  (9),  as 


dE 


d  (AP) 


For  P  =  2  atm. ,  the  maximum  error  for  excess  second  virial  coefficients 
is  0.06  cc/gm-mole  at  25°C  which  means  that  the  pressure  reading  error 
at  low  temperatures  is  relatively  lower  than  that  at  high  temperatures. 
This  will  give  the  maximum  error  of  3  percent. 

3.  The  largest  change  in  volume  of  the  oil  manometer  movement  is 
0.03  cc.  This  will  cause  the  pressure  change  of  0.6  mm  oil  and  it  may 
introduce  an  error  of  less  than  1  percent  in  the  pressure  reading. 

4.  The  error  due  to  the  impurity  in  the  pure  gas  is  negligibly 
small. 

5.  The  error  due  to  the  solubility  of  gases  in  the  manometer  oil  is 
negligible . 

6.  The  relative  error  in  reproducibility  of  the  experiment  is 
between  0.1  and  4  percent. 


APPENDIX  D 


DERIVATION  OF  EQUATION  USED 


D-l 


Derivation  of  Equation  (9) 


Write  the  pressure  of  a  gas  P,  in  terms  of  B,  C  and  density  d. 


P  =  RTd (1  +  Bd  +  cd  ) 


(D.l) 


Thus  for  the  ratio  of  final  to  initial  pressure. 


Pf  df (1  +  Bfdf  +  Cf df  ) 
d. (1  +  B.d.  +  C.d.2) 

l  li  li 


(D.2) 


since  pressure  in  the  bulbs  before  mixing  are  equal, 


P  =  dl(l  +  Bud  +  Culdl  ) 


d2('1  +  B22d2  +  C222d2  ^ 


(D.3) 


Neglecting  the  effect  of  the  pressure  change  on  the  density  of  the  gas 


that  is  at  room  temperature, 

then  d,  +  d. 


(D.A) 


The  mole  fractions  y,  and  y„  of  the  mixed  gases  may  be  expressed  in 

1  z 


terms  of  d^  and  d ^ ' 


1  dl  +  d2 


dl  +  d2 


D-2 


so , 


B, 


d1+d2 


B  +  2  -L-± -  B  + 

(d1+d2)z  12 


dl+d2 


B 


22 


substitute 


and 


since 


B 


12 


2  (B11  +  B22>  -  E 


then , 


B, 


f  d1 

d„d„ 

r  d„  i 

1 

1  2 

T5  _i_  0  t?  _L_ 

2 

di  d-d0 

11  J  n2  E  1 

d ,  +d  „ 

1  2 

v  J 

(d.+d^) 

[  1  u2  J 

B 


22 


(D.5) 


dl+d2 


C1U  +  3 


+  3 


r  > 

d. 

Z 

f  \ 

d. 

1 

2 

d  +d0 

d-+d0 

1  2j 

1  2 

^  / 

f  d,  ^ 

f  d 

z 

/•  "S 

d„ 

1 

2 

C122  + 

2 

d-  +d_ 

di  +d0 

d,+d0 

1  2i 

v.  ,7 

1  2j 

1  2j 

112 

3 


'222 


(D.6) 


AP  = 


Pf-P 


Thus 


AP 

P 


d.(l  +  Brdr  +  Crd  ) 
t  r  r  it 

V1  +  Bndi  +  cmV> 


-  1  (D  .  7) 


Substitute  Equation  (D.4),  (D.5),  and  (D.6)  into  (D.7)  and  after 


rearrangement , 


D-3 


AP 

P 


(d7~d  ) 

~u77 


d2E 


2(1  +  Budl  +  Cmd  ‘ 


where 


Clll(dld22  -  dl3)  +  C222(dl2d2  -  d23)  +  3dld2(dlGl+d2G2> 


8dl(1  +  Blldl  +  Cllldl  > 


G.  =  C.  .  .  -  -r  (2  C  .  .  +  C.  .  . ) 
1  1x3  3  m  333 


Rewrite  equation  (D.8) 


(D.8) 


E  =  2 


AP 

P 


<drV 

4dld2 


RTdld2 


(dl  -d2  ) (c111d1~c222d2) 

4dld2 


jHd^  +  d2G2)  (D  .9) 


substitute  for 


d. 

i 


P_ 

RT 


B.  . 

li 


RT 


+  (2B 


ii 


ciii> 


RT 


(D.10) 


Thus  equation  (D.9)  becomes 


(drd2>' 

dld2E  =  24PRT--^-3- 


/V  \ 


.  RT . 


+  4  ^dl  "d2  ^  (Cllldl  C222d2) 


-  |  did2(diGi  +  d2G2) 


(D.ll) 


D-4 


Then, 


(drd2) 


(B22"B11) 


'p 


RT 


+  (“2B22  +  C222  +  2B11  Clll^ 


rP 


,RT, 


+  ... 


(drd2>  = 


f' D 


iRT ; 


+  2  (B22-Bn)  ("2B22  +  C222 


+2B11  -Clli> 


RT 


+  ... 


(dl2-d22>  ’  2(B22-B11) 


^P  ^ 


iRT  i 


+  (5B112-5B222  -  2C111  +  2C222> 


rP  1* 


iRT  i 


+  ... 


dld2 


'p  12 


(B11  +  B22) 


rp  ,3 


RT , 

V.  ^ 


+  (2B112  +  2B222  +  B11B22 


Clll  "  C222) 


RT, 

V.  ^ 


+  (B11C222  +  B22cill  2B11B22 


"  2B11  B22> 


,RT, 


+  ... 


(<W 

dld2 


(B22  -  Bll> 


RT 


+  (B11  "  B22)(3B22  2C222  ”  3B11 


+  2C111) 


'P  13 


RT, 

v.  y 


+  •  .  . 


D-5 


dld2^dlGl  +  ^2G2^ 


(G1  +  G2) 


3 

p 

[rtJ 

G1(2B11  +  B22}  +  G2(B11  +  2B22} 


P_ 

RT 


+  ... 


(dl  d2  )(Cllldl  "  C222d2) 


2(B22  B11)('C111  "  C222) 


fr  ^ 


,RT; 


+  ... 


substitute  these  equations  into  equation  (D.ll) 


'p  '2 


RT . 

V  ^ 


'  (B11  +  B22> 


r p  i3 


RT 


+  (2Bn  +  2B22  +  BnB22  ”  CH1  “  C222) 


fP  ^ 


RT 


+  .  .  . 


2AP  _  1 
RT  2 


(B22  “  Bll)  RT 


+  (B11  _  B22> 


(3B222  -  2C222  -  3B112  +  “ill* 


fp 

3 

f  \ 

P 

IrtJ 

+  •  e  . 

RT 

2('B22  "  B11^C111  "  C222) 


rP 


RT 


+  ... 


r 


2 
4  * 


rp 


(G1  +  V 


'p  13 


tRTj 


G1(2B11  +  >22>  +  G2(B11  +  2B22) 


RT 


+  ... 


(D .  12) 


D-6 


rearrange  equation  (D.12) 


r 


2RTAP 


+  <BU  +  B22)  E 


E==+  1-|(B22-B11) 


"  4  (G1  +  g2) 


'p  1 


^RTj 


+  \ 


-  E(2BU2  +  2B222  +  B  Bm  -  Cnl  -  C222) 


2(B11  B22)(3B22  2C222  "  3B11  +  2Clli:) 


2(B22  "  Bll^  (C111  "  C222) 


+  4  g1(2B11  +  B22)  +  G2(B11  +  2B22) 


fP  13 


iRT ; 


Thus 


+ 


E  =  2RTAP  _  2 


(B11  ”  B22)  /4  "  E(B11  +  B22)/2  +  3(G1  +  G2)/8 


r  \ 

P 


RT 


”  2  {4(B11  B22^('C111  B11  )  (C222  B22  )  ^ 


2 


(2B11  +  B11B22  +  2B22  _  Clll  -  Clll) 


-  8  [G1(2B11  +  B22^  +  G2(B11  +  2B22>3, 


P  I2 


+  .  .  . 


(9) 


APPENDIX  E 

DETERMINATION  OF  SECOND  VIRIAL  COEFFICIENT  OF  CARBON  DIOXIDE 


E-l 


The  virial  equation  of  state  can  be  expressed  by  an  expansion 
using  the  pressure  as  the  independent  variable, 


PV 

RT 


1  +  B'P  +  C’P2  +  .  .  . 


using  the  first  two  terms  of  the  expansion 


PV 

RT 

1  +  B’P 

where 

B'  = 

B 

RT 

so 

PV 

1  +  2L 

RT 

1  +  RT 

or 

V 

S+B 

In  this  work,  volume  is  constant  and  equal  between  the  two  systems, 
then 


RT 

— —  +  B 

P,„  He 
lHe 


RT 

^2He  +  B2He 


(E.l) 


RT 

TZ  +  BlHe 

lHe 


RT1 

~C02  +  BlC°2 


(E.2) 


RT, 


RT. 


P2C02  +  B2C°2  =  P1C02  +  BlC°2 


(E .  3) 


P  =  P 
lHe  ICO, 


and 


(E .  4) 


E-2 


Equation  (E.l),  (E.2),  (E.3),  and  (E.4)  can  be  reduced  to 


RT, 


2He 


RT 

P7~  +  BlHe  "  B2He 
lHe 


>  =  P  +  AP 

2C0.  2He 


Thus 


RT. 


RT, 


B 


2C0, 


+  B 


1C0„  2C0, 


ICO, 


E-3 


TABLE  XXI 

Raw  Data  For  Second  Virial  Coefficient  of  Pure  Carbon  Dioxide 


T(°K)  P(mm  Hg)  AP(mxn  oil) 


Bco  (cc/gm-mole) 


Reference  point 

273.15 

1429 

263.15 

-21.47 

253.15 

-46.05 

-150.00 

-163.50 

-180.07 


Reference  point 

253.15 

1336 

243.15 

-26.50 

233.15 

-44.25 

-180.07 

-197.75 

-210.89 


Reference  point 

233.15 

1259.5 

223.15 

-33.35 

213.15 

-70.10 

-210.89 

-244.12 

-262.06 


Values  of  B„  Used: 
He 


T(°K) 

273.15 

253.15 


B  (cc/gm-mole) 


11.825 

11.880 


233.15 


11.950 


